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We derive the complete spectrum of gravitational waves induced by primordial scalar perturba- 
tions ranging over all observable wavelengths. This scalar-induced contribution can be computed 
directly from the observed scalar perturbations and general relativity and is, in this sense, indepen- 
dent of the cosmological model for generating the perturbations. The spectrum is scale-invariant on 
small scales, but has an interesting scale-dependence on large and intermediate scales, where scalar- 
induced gravitational waves do not redshift and are hence enhanced relative to the background 
density of the Universe. This contribution to the tensor spectrum is significantly different in form 
from the direct model-dependent primordial tensor spectrum and, although small in magnitude, it 
dominates the primordial signal for some cosmological models. We confirm our analytical results by 
direct numerical integration of the equations of motion. 
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1. INTRODUCTION 



Arguably the most striking prediction of inflationary cosmology [l| is the causal generation of nearly scale-invariant 
spectra of both scalar (energy density) and tensor (gravitational wave, GW) perturbations. The natural prediction is 
that the scalar and tensor amplitudes are comparable within one or two orders of magnitude of one another by virtue 
of the fact that both are created by the same de Sitter quantum process. The existence of a scalar spectrum is now 
firmly established by measurements of the cosmic microwave background (CMB) and large scale structure [3| , and 
its amplitude is well-determined. Tensor fluctuations, on the other hand, have yet to be detected, although current 
measurements have only begun to probe the expected range of amplitudes. 

Detecting primordial tensor fluctuations is an important milestone because it rules out a whole class of alterna- 
tive cosmological scenarios, like the ekpyrotic [4| and cyclic models [5(, which produce virtually the identical scalar 
spectrum as inflation but a completely different tensor spectrum. In particular, the primordial tensor contribution in 
ekpyrotic/cyclic models is exponentially smaller and more blue @. Detection of a primordial tensor signal is therefore 
widely regarded as a smoking gun signature of inflation. However, failing to detect the tensor modes at the expected 
level does not necessarily rule out inflation. The inflationary tensor signal can be suppressed by extra fine-tuning of 
the inflationary model and/or the addition of extra fields (e.g. hybrid inflation 0]) so that the background equation 
of state of the Universe, instead of changing smoothly during the final stages of inflation, undergoes a sequence of 
jerks and gyrations Q . A number of studies have discussed the limits to how far a search for the tensor spectrum can 
go based on detector sensitivity and foregrounds [9(. 




In(k) 



FIG. 1: Spectra of first- and second-order gravitational waves: This schematic illustrates the conjectured form of 
HgwW, the fraction of the critical density in gravitational waves per log-interval of wavenumber k, as derived in section [3] 
The topmost curve represents the typical first-order inflationary tensor spectrum. With fine-tuning, it can be suppressed below 
the level of the second-order, scalar-induced tensor perturbations (bottom curves). The bottom curves represent a sequence of 
times: matter-radiation equality (a cq ), redshift z — 100, and today (ao). The scalar-induced tensor spectra shown here are for 
a perfectly scale-invariant scalar input spectrum (n s — 1). If the scalar spectrum is blue (n s > 1) the induced tensor spectrum 
is enhanced on small scales (large k), while a red spectrum (n s < 1) suppresses tensor fluctuations on small scales (see section 
[5] for cautionary remarks about extrapolating spectra to very small scales using the large-scale power law form of the scalar 
spectrum). Ogw is of course ill-defined on superhorizon scales. On superhorizon scales (dashed lines) we therefore formally 
define the rescaled tensor power spectrum, k 2 Ph(k), but do not interpret it as an energy density of gravitational waves (see 
section [3J • 



At second order in perturbation theory the observed scalar spectrum sources the generation of secondary tensor 
modes [loj . In this paper, we analyze the stochastic spectrum of second-order gravitational waves induced by the 
first-order scalar perturbations. Since the scalar spectrum is already measured, this contribution to the tensor 
spectrum must exist and must be the same for both inflationary and ekpyrotic models because their predictions for 
the scalar spectrum match. For inflation, this second-order contribution is generically negligible, orders of magnitude 
smaller than the first-order contribution except for models with extreme fine-tuning. For ekpyrotic and cyclic models, 
the scalar-induced second-order contribution computed here is actually the dominant contribution on astrophysical 
and cosmological scales, because the first-order tensor spectrum is always exponentially small compared to the scalar 
spectrum. Hence, the calculation here supersedes previous predictions of the tensor spectrum for ekpyrotic and 
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cyclic models [6|. Because the gravitational wave spectrum we compute here is purely a consequence of the observed 
scalar spectrum and general relativistic evolution, any mechanism that accounts for the observed spectrum of scalar 
fluctuations also generates the same secondary tensor spectrum, provided Einstein's equations hold. Hence, this 
second-order signal provides an absolute lower limit on tensors from the early Universe. 

Our work builds on important earlier work by Mollerach, Harari and Matarrese [ill ] and Ananda, Clarkson and 
Wands [e|. Mollerach et al. 11] computed the effect of second-order gravitational waves on large-scale CMB 
polarization. They found that the second-order tensors dominate over the first-order signal if the ratio of the 
tensor-to-scalar amplitude on the largest observable wavelengths is r < 10 -6 . Then, more recently, Ananda, Clarkson 
and Wands [13] numerically studied the present spectrum of gravitational waves on very small scales accessible 
to direct detection experiments like the Big Bang Observer (BBO). Typically, the signal is expected to be at 
the extreme limit of the predicted sensitivity of BBO. Here we compute the complete spectrum of scalar-induced 
gravitational waves on all scales and discuss how it evolves with time. We analytically determine a critical scale 
above which second-order gravitational waves do not redshift. This non-trivial transfer function for second-order 
gravitational waves leads to an interesting feature in the current spectrum (see the schematic in Figure [IJ with 
a factor 10 7 enhancement of modes of order the horizon size at matter-radiation equality relative to the scale- 
invariant small scale spectrum. We confirm our analytical findings by numerical integration of the equations of motion. 

The outline of the paper is as follows: In we derive the evolution equations for second-order tensor fluctuations 
sourced by first-order scalar fluctuations. Allowing for an anisotropic stress contribution to the energy momentum 
tensor, we derive ageneral expression for the power spectrum of scalar-induced gravitational waves. This generalizes 
the work of Ref. [121 ]. In $31 we analyze the spectrum using various approximations and scaling arguments. These 
analytical estimates are confirmed by direct numerical integration of the equations of motion in <21 Finally, we 
discuss the implications of these results in $5] In two appendices we recall the Green's functions for gravitational 
waves and the transfer functions for first-order scalar fluctuations [l3j |. 

We use the following conventions: Throughout we employ natural units, H = c = 1, and (reduced) Planck mass 
Mp 2 = 8nG = n 2 , as well as 'East coast' signature for the metric, (—,+,+,+). Greek indices, fi, v = 0, . . . , 3, 
denote four-dimensional spacetime indices, while roman indices, i,j = 1, ... ,3, are reserved for spatial indices. The 
parameter rj is conformal time, a{rf) drj = dt. 



2. SECOND-ORDER TENSORS FROM FIRST-ORDER SCALARS 



Let us recall some basic facts about second-order perturbation theory, before deriving the explicit form of the evo- 
lution equations for second-order, scalar-induced tensors. We consider perturbations to a flat Friedmann-Robertson- 
Walker (FRW) background, g$ — a 2 (r/)r/^, 

9p.v = 9$ + $9»v , (1) 
where g$ satisfies the Oth-order Einstein equations, G$ = k 2 T$ , 

W 2 = ^V ( 0)) W 2_ W / = ^! (p (0) + p(0) ); U = drilna , (2) 

Here and P<°) are the homogeneous background density and pressure, respectively, and (...)' denotes a derivative 
with respect to conformal time, r\. Including to linear order the small quantum perturbations to the metric and energy 

density, the solution to the lst-order Einstein equations, G^ be decomposed into independent scalar, 

vector and tensor modes. At linear order, different fc-modes in Fourier space are independent. This is in contrast to 
the 2nd-order Einstein equations, G$ = k 2 T^) , where different fc-modes mix and scalar, vector and tensor modes are 
not independent. However, it is important to notice that, at second order, there is no mixing between second-order 

(2) 

scalar, vector, and tensor modes. On the other hand, there is a second-order contribution to the tensor mode, , 
that depends quadratically on the first-order scalar metric perturbation. This contribution, the "scalar-induced" 
tensor mode, is the focus of this paper. 1 



1 Scalar-induced vector modes were studied in |14| . 
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2.1. Evolution Equations 

To compute the second-order, scalar-induced tensor mode we begin with the following perturbed metric 



ds 2 = a 2 (r)) 



1 + 2$ (1) + 2$ (2) ) dri 2 + 2V7 2) dr?d 



1 - 2* (1) - 21< (2) ) 5ij + ^h i:j 1 dx'dx 1 



(3) 



(2) 

where hij = fiL and we have ignored first-order vector and tensor perturbations. Here and in the following, the super- 
scripts are formal labels for the order of the perturbation. The second-order Einstein tensor and energy-momentum 
tensor are [l5[ 



G 



(2)i 



and 



T 



+ (V 2) , $(2), T// 2 ) terms) + (diagonal part) <5j 

( 2 >*. = (>) + p(o>) w «v (1) + pWn^ . + p(D n (1) \. + p< 2 ^ , 



(4) 



(5) 



where p, P, w and II are energy density, pressure, velocity and anisotropic stress, respectively. We act on the spatial 
components of the Einstein equations with the projection tensor "Xy m [12j, 

1 ij ^Im ~ K 1 ij ± lm • 



(6) 



We will define the operator explicitly below, but we note here that it extracts the transverse, traceless part of 

any tensor and eliminates the terms involving $( 2 ), vpt 2 ', y^ 2 ', p( 2 ) and the scalar and vector parts of Yl^ 1 ■ in the 
second-order Einstein equations. Using the following first-order relations, 

P« EE C 2 P « 

2 



P 



(i) 



2 dA^'+n^A, 



K 2 a 2 (pW +p(o)) 



n 



1 



K 2 a 2 P(°) 

the evolution equation ([6]) can be written as follows, 



(7) 
(8) 
(9) 

(10) 
(11) 



(2)i 

where we have neglected the tensor part of IT • and defined 



Sij = 2$d i d j 3> - 2^6%® + 488*8 jib + d^dj® - 8^8^ - 9**6^$ + W^d^ 

Here, w = / p(°\ $ ee and '5 ee We define the Fourier transform of tensor metric perturbations as 

d 3 k 



%( x ,?/) = 



F e lk x [^k(»?)ey(k)+^k(»?)5tf(k)] 



(12) 



(13) 



(2tt) 3 / 2 

where the two time-independent polarization tensors and may be expressed in terms of orthonormal basis 
vectors e and e orthogonal to k, 

ey(k) ee -^=[e»(k)ej(k) - ej(k)ej(k)] : 



V2 l 

e«(k) = -^=[e 4 (k)e,(k)+e 4 (k)e,(k)]. 



(14) 
(15) 



5 



In terms of these polarization tensors, the projection tensor in ([6]) and (fTTj) is 

d 3 k 



7' lm c 
ij O/l 



„ik-x 



where 



(2tt) 3 / 2 



<Sz m (k) 



• x [e y -(k)e im (k) + e ti (k)e im (k)] 5, m (k) . 



e- 4kx 5 im (x'). 



(16) 



(17) 



(2tt) 3 / 2 

In Fourier space, the equation of motion for the gravitational wave amplitude h (for either polarization h or h) becomes 

h!i + 2Hh' k + k 2 h k = 5(k, ?7) , (18) 
where the source term, 5, is a convolution of two first-order scalar perturbations at different wavenumbers, 

S(M) = -4e Zm (k)S im (k) (19) 

= 4 



d 3 k 



(2tt) 3 / ; 



■e /m (k)fc ; fc„ 



7 + 3w 2c 2 



3(1 + iv) w 



2de 



*kfa)* k -kfa) + 11- ^ %(»>)*„_£(»?) 



2 C n i 



3(1 + ™) «i 



(20) 



Equation ([20]) reduces to the expression in [III in the limit * -> $, to -> 1/3 and c 2 -> 1/3. The limit * -> $, u> -> 
and c 2 — > was discussed in HHl. 



2.2. Power Spectrum 

The power spectrum of tensor metric perturbations, Ph(k,rf), is defined as follows 

(h k ( V )h K ( v )) = ^-6(k + K)P h (k, v ). 



(21) 



We now derive an expression for the power spectrum of second-order gravitational waves by solving equation (|18p . 
It is convenient to remove the Hubble damping term in (|18p by defining ah^ = «k, where t>k satisfies the following 
equation of motion 



k 



v k = aS . 



(22) 



The particular solution of ([18]) is then found by the Green's function method 

1 



a{rj) 



(23) 



where 



(24) 



Exact solutions to (|24|) for both matter and radiation domination are derived in Appendix [XI Substituting the solution 
(|23|1 into the expression for the tensor power spectrum (j21|) we find 



(hk(rj)hit(v)) = 



a 2 (rf) 



no 



d?72 / d?7ia(?7i)a(772)5k(??;7?i)gK(77;»?2) (S(k,fji)S(K,fj 2 )) 



(25) 
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The source term (l20l) may be written in the following form 



S(k,r)) = I d 3 ke(k,k)/(k,k^)V k _ £ ^, 



where 



and 



e(k,k) = e tf (fc)jfejjfej = fc 2 [l — a* 2 ] , /z 



kk 



(26) 



(27) 



/(k,k,7?) 



7 + 3w 2c 2 



.2c| 
w 



3(1 + u>) u> 

fc 2 N 



$(fer?)$(|k - icjr?) + 1 



21.2 



2c|fc_ 



*(fc7?)*(|k-k|77) 



1 



3H 2 



$(fctj)*(|k-k|T7) 
4 



8 



3(1 + w) m JH 
*'(ifc7?)*'(|k - k|r?) 



— ] i*(M*'(|k-k|?7) 



(28) 



Here we have split the first-order quantities into transfer functions, ^(krj), ^>(kr)), and primordial fluctuations ipk, 

$ u (?7) = <t>(k V ) Vk, *k(r?) = %{kri) fa . (29) 
The primordial fluctuations are characterized by the power spectrum, 



<V>kVk> = ^rfW(k + k) 



(30) 



k/ fc3 

Observationally, it is found that P{k) is nearly scale-invariant, so that the following parameterization is appropriate 

n B — l 

(31) 



P(fc) = -A 2 ,(fc ) - 



where recent CMB and large scale structure results @, Q imply A^(fc = 0.002 Mpc 1 ) = (2.40 ± 0.12) x 10 -9 and 
n s (k ) ~ 0.94 — 1.10 (0.95 ± 0.02; no tensors). Finally, the correlator in equation (f2"5)) can be computed using Wick's 
theorem 



(32) 



(5(k,%)5(K,^)) = / d J ke(k,k)/(k,k,^) / d d Ke(K,K)/(K ) K,%)(V' k _ £ V£V'K-KV'K/ 



6(k + K) / d 3 ke(k,k) 2 /(k,k,7y 1 ) /(k, k, t? 2 ) + /(k, k - k, fj 2 ) 



P(\k-k\)P(k) 
|k-k| 3 fc 3 



and the power spectrum of second-order gravitational waves is 



dk I dfi P(\k - k\)P(k) F(k,k, fi;ri) , 



(33) 



where 



[i- M 2 ] 2 fc 3 fc 3 n 



a 2 (v) Ik -k| 3 



dfj 2 dfj 1 a(fji)a(ri2)gk(V]Vi)9k(mV2)f(^,K'ni) /(k k fj 2 ) + f(k, k - k, fj 2 ) ■ (34) 



Notice that the power spectrum Ph(k,if) is defined completely in terms of the Green's function (Appendix [X]), the 
transfer functions $ and \& (Appendix [B]), and the primordial power spectrum of first-order scalar fluctuations, P(k) 
(WMAP 0). 
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3. ANALYTICAL DESCRIPTION OF THE SPECTRUM 



In this section we estimate the complete spectrum of scalar-induced gravitational waves analytically. To simplify 
the analysis we neglect anisotropic stress and set ^ = $. In section Q] we evaluate the exact spectrum numerically 
including anisotropic stress and show that this gives only a small correction. With \& = <J>, the source term of the 
equation of motion (TT8|) can be expressed solely by the Bardeen potential <f>, 



h'l + 2Hh' k + k 2 h k = S($(kr))) , (35) 
and /(k, k,??) in equation (f28|) is expressed by a single transfer function 

3 4 f(K k, rj) = 2(5 + 3u>)$(|k - k|77)$(|k|77) + 4 (2r/$(|k - k|jj) + ?? 2 $'(|k - k|»j)) S'flkfo) . (36) 
In Appendix [B] we show that the transfer function for first-order scalar modes can be written in the following form 

*(*»?) =\^ltT w 

l +fe 2 t? 2 V > Veq 



Amplitude 



S c /k 2 - 



h D =0 



k h (i) = S (i) /k 2 




S<% 2 = h (,) 



a k 



ln(a) 



FIG. 2: Evolution of scalar source and induced gravitational waves. Second-order tensors, h, are generated when the 
mode k enters the horizon at au ■ If horizon entry occurs during the radiation dominated era, then the scalar source decays as 
a -7 until matter-radiation equality, a oq . During matter domination the scalar source terms remains at a constant value, S^ f . 
Gravitational waves redshift like a -1 as long as h > 5' f '/fc 2 , but remain at a constant amplitude maintained by the constant 
source term after that, a > a%. 



To study the generation of h induced by S we make the approximation that gravitational waves are produced 
instantaneously when the relevant mode enters the horizon. The subsequent evolution of the tensor mode is scale- 
dependent and determined by the time evolution of the scalar source term (see Figure [2]). Scalar- induced gravitational 
waves redshift as long as their magnitude is greater than S/k 2 . After that they freeze at a constant value maintained 
by the constant source term during matter domination. We define the transfer function for scalar-induced gravitational 
waves, t(k,rj), as follows 

hkirfi^tik^hf , (38) 

where is the value of just after the instantaneous generation of gravitational waves after horizon entry (see 

Figure [J). We estimate by dropping time derivatives in the equation of motion (|35[) (since kr\ > 1 after horizon 
entry) 



4° - ^ 



(39) 
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In £|3- II we calculate the initial power spectrum at the time of horizon crossing, 

PjtHkMty^^siihP) 2 ), (40) 

where r)i(k) ~ fc -1 is the conformal time when a comoving scale fc enters the horizon. This initial spectrum is processed 
using the tensor transfer function, t(k,rf), which we derive in £13.21 Finally, in £)3.31 we put these results together and 
compute the relative energy density of scalar-induced gravitational waves 



2 cqk cq 



u2 -t\k,r,)P^(k). (41) 



3.1. Power Spectrum at Horizon Crossing 

In this section we estimate the fc-scaling of the horizon power spectrum of scalar-induced gravitational waves. The 
horizon amplitude is estimated from the equation of motion (|35p as follows 

h k ~ ^« ~ iy d 8 kfc 2 (l-^)$(fc»»)$(|k-k|»»)^_£, (42) 

and its power spectrum is, 

~ p^a |d 3 kd 3 Kp£ 2 (l- M 2 )(l-A 2 )<^^ 

~ ^(k + K) / d 3 kfc*(l - ^)H 2 (kVi)* 2 (\*- khi)^ P f ~- k 3 D ■ (43) 
fc 4 J fc 3 |k-k| 3 

Hence, 

if (fc) ee |^<(4 i) ) 2 > ~ J / dW(l - M 2 ) 2 $ 2 (^)$ 2 (|k - kM^ ^ k ~-jf • (44) 

To compute (|44[) we use the transfer function for the scalar potential (|3T[) and assume a scale-invariant spectrum, 
P(fc) = fA 2 ^. (The scale-dependence of the scalar spectrum can be reinserted at the end of the computation). Here 
we have defined A 2 ^ « 10 -9 as a measure of scalar power on COBE scales fco ~ 0.002 Mpc -1 . Hence, for the radiation 
dominated phase, we have, 

PP&Mk)) ~ / dfc / d^l-/^ 2 - l - : 



fc Jo J-l (fc 2 + fc 2 - 2fcfc^) 3 / 2 [l+(fc/fc) 2 ] 2 [l + (fc 2 +fc 2 -2fcfc/x)/fc 2 ] 2 

X 3 1 1 



/•OO /■ 1 

A 4 , J dx J ^dMl-M 2 ) 2 - 



(1 + x 2 - 2xfi) 3 / 2 (1 + x 2 ) 2 (2 + .t 2 - 22^) 2 
cx A 4 ,, (45) 

where we defined x = &/&;. On the other hand, for the matter dominated phase, we have, 

^^17 f -\7 f 9\9 1 1 



p£\kMk)) ~ dkj Mi-^f 



(fc 2 + fc 2 - 2fcfc/x) 3 / 2 [1 + (fc/fceq) 2 ] 2 [1 + (fc 2 + fc 2 - 2kkn)/k, 



2 12 



f°° /- 1 a; 3 1 1 

= A* / dcy^ d/i(l - /i 2 ) 2 (i - ^ _ — yj- 2 (1 - ^ y2)2 [1 + {1+x2 _ 2xfl ) y 2]2 ' ( 46 ) 
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where x = k/k and y = k/k eq <§; 1. Neglecting the /i-terms we find 

^3 



dx- 



x 



1 



(l + x 2 ) 3 / 2 (1 + xV) 4 
5(1 + 6y 2 )arccosy 81 + 28y 2 - 4y 41 



16y(l - y 2 ) 9 / 2 



48(1 -y 2 ) 4 



5^ 
A 4 ^£a 

The power spectrum at horizon crossing therefore scales as follows 



(47) 



P«(fc)ocA^ ^ 



(48) 



3.2. Second-Order Tensor Transfer Function 



To compute the transfer function for second-order, scalar-induced gravitational waves we need to estimate the time 
evolution of the source term. In particular, the transfer function for modes that enter the horizon during radiation 
domination is sensitive only to ratio of the source terms at horizon crossing, <SW and the asymptotic value after 
equality ( see Figured]). We parameterize the decay of the source term during radiation domination as follows 



Ctk 



7(fc) 



(49) 



where we have allowed for a scale-dependence of the effective decay rate. In the following we put limits on j(k) by 
considering the asymptotic evolution of subhorizon modes (krj 3> 1). 



Let x = |k — k| 77 and y = |k|^. Using the following relations 



\ krj J 



= 1 



V krj 



V krj 



and 



the source term may be written as follows 



5(k,7?) oc — / dy d/x [1 
V Jo J-i 



,,5/2 



x(y,/i) 3 / 2 



f{x,y) ■ 



where 



2f{x,y) = $(x)$(y) 6 - y 2 $(x)$(y) 



6 



ly 2 



1 



ly 2 



1 



(50) 

(51) 

(52) 
(53) 



Let us estimate the integral (f5Tl) . The limit y — > is clearly suppressed by the phase space factor y 5 / 2 in the integrand. 
The limit x — > (y — > fery, — >• 1) is suppressed by the projection factor [1 — fi 2 ]. To see this, first take the limit 

y -> 



f- 

V krj 



2[l-/x] 



1-M 2 

3-3/2 



(1- M )V4 (1+M) . 



(54) 



This shows that the integrand vanishes in the limit x —* 0. In addition, large x and y are suppressed by the transfer 
function f(x,y) (i.e. the decay of the Bardeen potential on subhorizon scales). The dominant contribution to the 
integral (f!)l"j) therefore comes from regions of phase space where k ~ k (y ~ fcry) and |k — k| ~ k (x ~ zx ~ 0). Let 
us therefore write 



Soc -i y dlny y dln(l- M )(l-/i) 2 (l + /x) 



7/2 1 



1 



2-3/2 x 2 _|_ 1 y 2 + } 



6- 



1 



?/ 



ly 2 



(55) 



10 



and take the subhorizon limit y, x — > kr] > 1, fi - 

c 1 (kr]) 7 ' 2 1 



r) 2 (/cr?) 3 / 2 [(fcr?) 2 + l] 2 



6 _ i k vf 



[{kri) 2 + lp 



' ' oc^. (56) 



2 (Wl2 ~ „4 



?y 2 (kr]) 



Hence the source term decays at most as a~ 4 after the mode k enters the horizon during the radiation dominated 
era, i.e. "f(k) < 4. In fact, we expect the source term to decay considerably slower than that for a while after horizon 
crossing. The source will decay more and more quickly as the horizon grows much larger than the wavelength of the 
mode, finally reaching the asymptotic behavior that is proportional to a~ 4 . This leads us to expect that the effective 
7 in equation (|49|) will be significantly smaller than 4. Numerically, we find 7 3 (see <j4|. 

The transfer function for second-order gravitational waves is considerably different from the transfer function for 
first-order gravitational waves. First of all, modes which enter the horizon during matter domination have constant 
source terms and hence do not decay 

t(k,ri) = l, k<k cq . (57) 

Next, consider the evolution of the scalar source term and induced gravitational waves for modes that enter the 
horizon during the radiation dominated era (see Figure [2]) ■ Here, k = akH defines the time of horizon entry (a*;) for 
a mode of wavenumber k. We assume that h grows very rapidly after horizon entry to become of order the source 
term. Then S decays as a -7 (where our previous discussion implies 7 < 4) while h redshifts as a -1 until h is equal 
to the final source term during matter domination at a* k . For a > ajjl, h stays constant. We therefore have, 

h® ~ a* k 5« ~ \a e J ' (&S > 

and find 

^ = (2hV = (JLY\ (59) 

a k V «oq^ \k oq S 

For a fixed time 77, subhorizon modes with sufficiently large k have never settled down. The critical wave number at 
a time tj can be obtained by substituting a* k = a(rj) into equation (|59p . 

W = (^) 1/(7 ~V (60) 



Modes with k > k c (r]) simply redshift like a , 

t(k, V ) = ^- = ^-^-, k>kM- (61) 

a (77) a(r)) kr] cq 

The transfer function for second-order gravitational waves therefore takes the following interesting form 

{1 k < k eq 

(£) 7 fc c q < k < k c ( v ) (62) 



3.3. Spectrum of Scalar-Induced Gravitational Waves 



Substituting the power spectrum at horizon crossing (|48[) and the tensor transfer function (|62[) into equation ([IT 
for the relative spectral energy density of gravitational waves at time 77, we find 



^ { Gw( k ^) 



A {2) A 



k \ 



( a(ri) k 

^ca ^cq 



k < k, 



cq 



a(v) ( k \ 



-(27-2) 



a (v) 



k cq < k < k c (rj) 
k > k c (rj) 



(63) 
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where an overall normalization constant A^, hasn't been fixed by our analytical arguments. In Sj4]we find AqL, w 10 
and 7 rs 3 (This is consistent with the normalization of the small-scale spectrum in 12]). Figure [1] summarizes this 
conjectured form of the scalar-induced gravitational wave spectrum. If the scalar spectrum can be treated by a power 
law with constant spectral index over a large range of scales, then a blue scalar spectrum (n s > 1) enhances the tensor 
spectrum on small scales, while a red spectrum (n s < 1) suppresses secondary tensor fluctuations. For comparison, 
the first-order spectrum of primordial gravitational waves in inflationary models can be expressed as 

"gw 



(M)=4VoA^o)(fr l f^fe) (64) 

where A^^ = 4.2 x 1CP 2 and ro = P p^-j is the first-order tensor-to-scalar ratio evaluated on the scale of today's 

horizon, k — ko w 0.002 Mpc -1 . For single-field inflation, the spectral index of the primordial tensor spectrum, n t , 
is related to the tensor-to-scalar ratio by the slow-roll consistency relation, n t = — ro/8. Matter-radiation equality is 
normalized by recent observations Q a e q ~ ao/3400. 

The current tensor spectrum (77 — t]q) on large scales (k = k e q) and on very small scales (k ^> fc cq ) (assuming 
inflation and n t « and n s ~ 1) satisfies 

r o — W\aT7iT ~ T7T r o (65) 



ng> = A ( H^n(ko) 10 



'GW^ ~~ "GW 

and 

"lf (fc>>fccq) « r /°f 1 - ~ 10 6 r . (66) 

Hence, for inflation, on large scales the second-order, scalar-induced contribution today, in fact, dominates over the 
first-order contribution. This reflects the fact that second-order gravitational waves do not redshift on large scales, 
while first-order gravitational waves redshift on all scales. On small scales the first-order contribution dominates 
unless ro < 10~ 6 . For ekpyrotic/cyclic models, the first-order contribution (due to direct quantum fluctuations of the 
metric) is suppressed at k — ko by 60 orders of magnitude compared to the inflationary signal and the spectrum is 
blue. Hence, in these models, the scalar-induced tensor modes, comprise the dominant contribution on all scales. 

CMB observations probe the time of photon-baryon decoupling at ocmb ~ 3 a oq and scales with k < fccMB = 

M 2 . . . 

fc eq . The gravitational wave spectrum at that time satisfies 



\ QCMB / 







gw(^ = ^-cmb) ^ -^gw v° CMB / ( fc cq \ _ 10 6 ro (67) 



and 



"|f (fc>>fccq) « r A™ . 2 ),. ~ 10 6 r ■ (68) 
^(fc»fc cq ) 4w A ^ (fco) 

Hence, at decoupling the first-order tensor signal dominates over the second-order signal if ro > 10~ 6 . (This is 
consistent with the result of Mollerach et al. who claim that second-order gravitational waves only have a 

significant imprint on the CMB if ro < 10~ 8 . This corresponds to second-order gravitational waves dominating 
over first-order gravitational waves at the time of recombination. Second-order effects can become visible for larger 
ro < 10~ 6 if late time polarization generated by reionization is considered [111].) 

For completeness, let us consider the power spectrum on superhorizon scales, e.g. during the matter dominated 
phase. On superhorizon scales, k <C 7i, the initial amplitude can be estimated from the equation of motion (|35p by 
ignoring the gradient term and approximating the time derivatives by factor of 7i, 

n 2 h^ <n ^S, (69) 
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while k 2 h^ ~ S estimates the initial amplitude on the horizon scale. Hence, the initial power spectrum, Pu\ on 
superhorizon scales is simply (k/Ti.) 4 times the spectrum on the horizon scale 

P h\k) OC (J£\ X A^(fco)^ CX fc 3 , k < fchor < fceq . (70) 

(2) 

Although the tensor power spectrum, Ph, is a well-defined gauge-invariant object on super- horizon scales, f^yy is not. 

In particular, equation (|4ip is only defined on subhorizon scales. Nevertheless, we formally define £1 ^ oc k 2 P^ ] oc fc 5 
on superhorizon scales, but do not attribute physical meaning to it. This definition is useful, since all our results 

(2) 

are presented in terms of ^q4? an d the shape of the superhorizon spectrum gives a simple consistency check for the 
numerical analysis. 



4. NUMERICAL RESULTS FOR THE EXACT SPECTRUM 



The spectrum of scalar-induced gravitational waves that we derived in $2] and discussed analytically in $3] can be 
evaluated exactly using standard numerical methods. The time evolution of the first-order perturbation variables 
necessary to compute the spectrum, p(ky) and ^(kn), is obtained from publicly available Einstein-Boltzmann codes 
such as CMBFAST [lj| or CAMB 17]. We first store the time evolution of $ and * in fc-space, then convolve 
them according to equation (|52")l . In practice, the range of k is taken be [10 -5 Mpc -1 , 500 Mpc -1 ] and variables are 
evaluated at 50 uniformly spaced points per log-interval of k. We have checked that our results are stable under 
variations of the fc-space boundaries and the discretization. 
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FIG. 3: Numerical spectra of scalar-induced gravitational waves (lower curves) and the scale-invariant primordial tensor 
spectrum for an inflationary model with tensor-to-scalar ratio r = 0.1 (upper curve). The scalar-induced spectra are shown 
at three different epochs, z + 1 = 3400, 100, and 1. Each curve has been extended, for pedagogical reasons, to modes with 
small wavenumbers k that lie outside the horizon at the given epoch (dotted range of the three lower curves). Note that 
current (z + 1 = 1) scalar-induced contributions cross the primordial inflationary contribution at intermediate wavelengths, as 
suggested by the schematic in Figure [T] The simulation assumes a flat ACDM cosmology with the following model parameters: 
A^(fc = 0.002 Mpc" 1 ) = 2.4 x 10~ 9 , n s = 1, n t = 0, r = 0.1, Sl b h 2 = 0.022, fl m h 2 = 0.11, h = 0.7. 



In the numerical analysis it is possible to incorporate the difference between $ and $ resulting from anisotropic 
stress of the fluid because the Boltzmann equations of photons and neutrinos are solved explicitly in the code by 
expanding their distribution functions into multipole moments. (Neglecting anisotropic stress from neutrinos implies 
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~10% errors for both first-order scalar and tensor perturbations [H, Ell; see Figure [5] in Appendix iBl For second- 
order tensors we find that the inclusion of anisotropic stress typically has less than 1% effect on the amplitude of 
the spectrum.) Finally, we should mention here that our definition of c 2 , equation (|7|), relates first-order pressure 
and energy density perturbations in the total matter, and p^ 1 ', including entropy perturbations, and thus it is 
scale-dependent. We have derived and incorporated this numerically. 

The numerically calculated spectra are shown in Figure [3l The shape of the scalar-induced gravitational wave 
spectrum agrees well with the analytical results of the previous section and the schematic diagram in Figure [TJ 
except for the fine-scale oscillations and the modest smoothing of the enhanced feature at large wavelengths. The 
scalar-induced spectrum is derived directly from observations of the scalar perturbations plus general relativity and 
is, in this sense, independent of the cosmological model for generating the primordial perturbations, e.g., inflation vs. 
ekpyrotic/cyclic. However, the transfer function does depend weakly on the expansion rate and composition of the 
universe, and, hence, the cosmological background parameters must be measured or otherwise specified. 



5. DISCUSSION 

Precise cosmological observations 0, Q have confirmed the existence of a nearly scale-invariant spectrum of pri- 
mordial scalar fluctuations. These scalar fluctuations induce a second-order contribution to the spectrum of tensor 
perturbations that must be present for any cosmological model that accounts for the observed scalar spectrum. In 
particular, the computation of the second-order gravitational wave signal does not assume that the primordial per- 
turbation spectra were generated by inflation - it only relies on the observed spectrum of scalar perturbations and 
general relativity. 

In this paper we have computed the scalar-induced spectrum of gravitational waves produced in the early Universe 
and used the Einstein equations to evolve it to the present. We have extended previous approaches to this problem 
[HI by considering the complete cosmic history for the evolution of second-order gravitational waves and allowing 
for anisotropic stress in our numerical work. 

Perhaps the most interesting theoretical feature is that second-order gravitational waves do not redshift on large and 
intermediate scales, but are maintained at a constant amplitude by the scalar source terms. This leads to a transfer 
function for the scalar-induced gravitational waves that produces a (nearly) scale-invariant spectrum on small scales 
and interesting scale-dependence on large and intermediate scales. In particular, there is a peak in the current 
spectrum of second-order gravitational waves at the scale of the comoving horizon at matter-radiation equality (see 
Figures [TJ [3] and [J} that is likely to exceed the primordial tensor spectrum at the present epoch. Unfortunately, there 
are no known methods for directly probing the present gravitational wave spectrum on these scales (corresponding to 
the size of superclusters today). At earlier times, such as recombination, the feature was much smaller and so only 
has small effects, e.g. on the CMB [ll[. Hence, this substantial feature is likely to remain of purely academic interest 
in the foreseeable future. 

On much smaller scales, which may be accessible to space-based laser-interferometer experiments (for a nice 
discussion see Ref. [2l|V there are no measurements of the scalar perturbation spectrum, so one must rely on 
extrapolating from what is known about scalar perturbations on large scales, e.g. from measurements of the CMB 
and large scale structure. Since these two wavelength regimes are separated by 16 orders of magnitude, extrapolation 
uncertainties can have important effects. For example, the dashed line in Figure^ which illustrates the extrapolation 
of the tensor spectrum based on a perfectly scale-invariant (n s = 1) scalar spectrum, is an estimate of the induced 
small-scale tensor signal. Assuming that a nearly constant spectral index is a valid approximation from CMB scales 
to the smallest scales, a blue spectrum (n s > 1) enhances scalar-induced tensor modes on direct-detection scales, 
while a red spectrum (n s < 1) suppresses them. One might imagine improving on this simplistic extrapolation 
by including running of the scalar spectral index, a s = dn s /d]nk, although observationally a s is only weakly 
constrained. However, there are interesting cases where this approach is inadequate [2l[ and can even lead to 
qualitatively misleading conclusions [22^ . Alternatively, for explicit inflationary models the whole spectrum can be 
computed directly from the inflaton potential V(<fi) without expanding with respect to the CMB scale. 2 In particular, 



where p, > Mp and <f> < p,, the exact computation 



for hybrid inflation models of the form V(<f>) — Vq 1 + (j^j 

predicts enhancement of the scalar-induced tensors on direct-detection scales, while the power law approximation 
(for p > 3) with n s > 1 and a s < predicts an overall suppression [22| . In a follow-up paper, we explore examples of 
this type that may produce large signals at BBO scales [241 ]. 



2 For the first-order inflationary tensor spectrum this approach was followed by |2:J . 
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FIG. 4: Observational prospects. Shown are the theoretical predictions for the present spectrum of primordial inflationary 
gravitational waves and scalar-induced gravitational waves, as well as current (solid bars) and future (dashed bars) experimental 
bounds (figure modified from [1]). The range of amplitudes for the primordial tensors corresponds to minimally tuned models 
of inflation as described in With more fine-tuning this signal can be suppressed. In contrast, the amplitude of the scalar- 
induced tensors is fixed by the observed amplitude of scalar fluctuations and therefore provides an absolute lower limit on 
the stochastic gravitational wave background. The CMB constraints depend on assumptions about the transfer function for 
gravitational waves to extrapolate constraints obtained at decoupling to the current spectrum. Since second-order gravitational 
waves do not redshift on CMB scales, the CMB observations imply separate constraints on the current first- and second- 
order spectra. These constraints are labeled CMB(l) and CMB(2), respectively. The dashed section of the scalar-induced 
tensor spectrum illustrates extrapolation from CMB to direct-detection scales using a scale-invariant scalar spectrum (n s = 1) • 
Important uncertainties in the extrapolation between CMB and BBO scales are discussed in the main text. 



Finally, let us consider the possible scenarios for future observations of the first- and second-order tensor signals 
and what they would signify: 

• If a spectrum of gravitational waves is observed that conforms to a nearly scale- invariant, first-order tensor 
signal and with r > 10~ 2 (as shown in Figure [4]), this would be a spectacular confirmation of the inflationary 
model of the universe and completely rule out ekpyrotic/cyclic models. The observation of such a signal with 
r < 10~ 2 would also rule out ekpyrotic/cyclic models but, in addition, inflationary models would be pushed into 
a fine-tuned region of parameter space Q. 

• If no nearly scale-invariant, first-order tensor spectrum is detected but the scalar-induced, second-order tensor 
spectrum is observed (either by extremely sensitive CMB polarization experiments or small scale direct detection 
interferometers like BBO), then inflation could only be compatible with extreme enough fine-tuning to suppress 
the first-order contribution to the tensor signal, and alternatives like the ekpyrotic/cyclic models would be 
favored. 

• If future experiments show that there is no tensor signal at or above the level of the predicted scalar-induced 
tensor spectrum, either general relativity or the interpretation of the scalar fluctuations would have to be amiss. 

In practice, observing the second-order tensor signal is a long way off, at best, and perhaps even impossible given 
our current understanding of astrophysical foregrounds and detector limitations for both CMB and direct detection 
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experiments [3, l20(. Nevertheless, we consider it interesting that the observed level of scalar fluctuations implies a 
model-independent lower limit on gravitational waves from the early Universe whose detailed features can be computed 
from a general relativistic description of cosmic evolution. 
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APPENDIX A: GREEN'S FUNCTION FOR GRAVITATIONAL WAVES 



In this section we derive exact Green's functions of equation (|18|) for both the radiation and matter dominated eras. 
The Green's function of a general second-order differential equation, Cg = S(i] — fj) is defined as follows 



vi{ti)v2{v) - vi(fj)v 2 (rj) 



Vi(fj)v 2 (fj) - vi(fj)v' 2 (fj) ' 



(Al) 



in terms of the two homogeneous solutions v\ and v 2 , which satisfy Cvi — 0, for a general differential operator C 
During the radiation dominated era the Green's function for the gravitational wave problem (|24|) reduces to 



9k + k2 9k = 8(r}-fj), 

which has the following homogeneous solutions 

V\ = sin(fc?7) , v 2 = cos(fc?y) . 
Hence, the Green's function during the radiation dominated era is 

9k(v, ff) — T[ s ^ n Q tr l) cos(fc^) — sin(fc?7) cos(fcTy)] 
k ' ' 



V < Vcq ■ 



During matter domination equation (|24p reduces to 



// 

9k 



which has the following homogeneous solutions 

vi=r)ji(x), v 2 = rjy-L(x), x = krj, 



(A2) 
(A3) 

(A4) 

(A5) 
(A6) 



where ji(x) and yi(x) are spherical Bessel functions. The Green's function during the matter dominated era therefore 
is 



9k(mv) 



XT 



j 1 (x)y 1 (x) - ji(x)yi(x) 



(A7) 



APPENDIX B: TRANSFER FUNCTION FOR FIRST-ORDER SCALAR MODES 

The first-order scalar perturbations $ and $ in equation (J3J) satisfy the following constraint equation [l3| 

fc 2 ($ - *) = -An 2 a 2 [p 7 9 2 + p v M 2 \ , (Bl) 

where Q 2 and N 2 characterize the quadrupole moments of the photon (7) and neutrino (v) anisotropics, respectively. 
02 and A/a are determine by the solution to the Einstein-Boltzmann equations. In practice, these are solved 
numerically using CMBFAST or CAMB [r| (see Figure [SJ. 

Since, 02 and J\f 2 are typically negligibly small, analytical studies often assume $ « ^. In this case, the first-order 
equation of motion for the Bardeen potential is (e.g. [l3j) 

+ 3(1 + c 2 s )H& - c 2 A<i> + (2H' + (1 + 3c 2 )H 2 ) cx SS , (B2) 



where the RHS is non-zero only in the presence of entropy perturbations, SS. For p = wp and in the absence of 
entropy perturbations (i.e. SS = 0) this becomes 

1 + 6w 77 

Equation (|B3[) has the following exact solution 

^(V) = y- a \Ci(k)J a (y) + C 2 (k)Y a (y)} , y = ^Ekr, , a =l(^±^L), (B4) 
L J 2 \ 1 + 6w J 
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FIG. 5: Transfer functions for <3> and fy. Neutrino anisotropic stress leads to 0(10)% difference between <!> and ^ [l£ 



where J Q and Y" a are Bessel functions of order a. During the matter dominated era (w = 0) this becomes 

$ k (,) = C 1 (fc) + ^P, 



whereas during the radiation dominated era (w = i) we find 



Ci(fc)( 



3/ 

sin y 



cos y 



C 2 (k)( 



cos y 
v V 



s\ny 



At early times y = yjwkr] <C 1 this becomes asymptotically the primordial value, 

lim $ k (??) = Ci(fc) = ^k, 



(B5) 



(B6) 



(B7) 



where we have dropped the decaying mode (C2 = 0) . For the growing mode solution we therefore obtain the following 
transfer function 



[ const. r\ > rj cq 

From this we see that superhorizon modes (krj -C 1) freeze during the radiation era 

i>(fcry) = 1 + O^knf) , kr] < 1 , 77 < ?y oq , 
while subhorizon modes (kr) > 1) oscillate and decay as a~ 2 



&(kn) = - 



cos[krj] fr/ k 



(kr)) 2 



cosffcr;] = — (— ) cos[fc?7] , krj > 1 , i] < r] cq . 
Tj / V a ' " '' 



(B8) 



(B9) 



(B10) 



Ignoring oscillations we therefore may write the following expression valid for both superhorizon and subhorizon modes 

1 



1 + {krif 

The Bardeen potential freezes on all scales during matter domination. 



(Bll) 



18 



[1] A. H. Guth, "The Inflationary Universe: A Possible Solution To The Horizon And Flatness Problems," Phys. Rev. D 23, 
347 (1981); 

A. D. Linde, "A New Inflationary Universe Scenario: A Possible Solution Of The Horizon, Flatness, Homogeneity, 
Isotropy And Primordial Monopole Problems," Phys. Lett. B 108, 389 (1982); 

A. Albrecht and P. J. Steinhardt, "Cosmology For Grand Unified Theories With Radiatively Induced Symmetry 
Breaking," Phys. Rev. Lett. 48, 1220 (1982). 
[2] D. N. Spergel et al. [WMAP Collaboration], "First Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: 
Determination of Cosmological Parameters," Astrophys. J. Suppl. 148, 175 (2003) [arXiv: astro-ph/0302209 ; 
H. V. Peiris et al., "First year Wilkinson Microwave Anisotropy Probe (WMAP) observations: Implications for inflation," 
Astrophys. J. Suppl. 148, 213 (2003) [arXiv:astro-ph/0302225] ; 

D. N. Spergel et al., "Wilkinson Microwave Anisotropy Probe (WMAP) three year results: Implications for cosmology," 
|arXiv:astro-ph/0603449| . 

[3] M. Tegmark et al., "Cosmological Constraints from the SDSS Luminous Red Galaxies," arXiv:astro-ph/0608632 

[4] J. Khoury, B. A. Ovrut, P. J. Steinhardt and N. Turok, "The ekpyrotic universe: Colliding branes and the origin of the 

hot big bang," Phys. Rev. D 64, 123522 (2001) [arXiv:hep-th/0103239] . 
[5] P. J. Steinhardt and N. Turok, "A cyclic model of the universe," Science 296, 1436 (2002); 

P. J. Steinhardt and N. Turok, "Cosmic evolution in a cyclic universe," Phys. Rev. D 65, 126003 (2002) 

[arXiv:hep-th/01 11098] . 

[6] L. A. Boyle, P. J. Steinhardt and N. Turok, "The cosmic gravitational wave background in a cyclic universe," Phys. Rev. 

D 69, 127302 (2004) [arXiv:hep-th/0307170] . 
[7] A. D. Linde, "Hybrid inflation," Phys. Rev. D 49, 748 (1994) [arXiv:astro-ph /9307002] . 

[8] L. A. Boyle, P. J. Steinhardt and N. Turok, "Inflationary predictions reconsidered," Phys. Rev. Lett. 96, 111301 (2006) 
|arXiv:astro-ph /0507455 . 

[9] J. Bock et al., "Task Force on Cosmic Microwave Background Research," arXiv:astro- ph/06 0lToTj 
[10] S. Matarrese, O. Pantano and D. Saez, "A General relativistic approach to the nonlinear evolution of collisionless 
matter," Phys. Rev. D 47, 1311 (1993); 

S. Matarrese, O. Pantano and D. Saez, "General relativistic dynamics of irrotational dust: Cosmological implications," 
Phys. Rev. Lett. 72, 320 (1994) [arXiv:astro-ph/9310036] ; 

S. Matarrese, S. Mollerach and M. Bruni, "Second-order perturbations of the Einstein-de Sitter universe," Phys. Rev. D 
58, 043504 (1998) [arXiv:astro-ph/9707278] ; 

H. Noh and J. c. Hwang, "Second-order perturbations of the Friedmann world model," Phys. Rev. D 69, 104011 (2004); 
C. Carbone and S. Matarrese, "A unified treatment of cosmological perturbations from super-horizon to small scales," 
Phys. Rev. D 71, 043508 (2005) [arXiv:astro-ph/040761l| ; 

K. Nakamura, "Second-order gauge invariant cosmological perturbation theory: Einstein equations in terms of gauge 



invariant variables," Prog. Theor. Phys. 117, 17 (2007) arXiv:gr-qc/0605108 



[11] S. Mollerach, D. Harari and S. Matarrese, "CMB polarization from secondary vector and tensor modes," Phys. Rev. D 

69, 063002 (2004) arXiv:astro-ph/031071l] . 
[12] K. N. Ananda, C. Clarkson and D. Wands, "The cosmological gravitational wave background from primordial density 

perturbations," |arXi v:gr-qc/0612013 
[13] V. Mukhanov, Physical foundations of cosmology, (Cambridge Univ. Press, Cambridge, 2006); 

S. Dodelson, Modern cosmology, (Academic Press, San Diego, 2003) . 
[14] S. Matarrese, S. Mollerach, A. Notari and A. Riotto, "Large-scale magnetic fields from density perturbations," Phys. 

Rev. D 71, 043502 (2005) [arXiv:astro-ph/0410687] ; 

R. Gopal and S. Sethi, "Generation of Magnetic Field in the Pre-recombination Era," Mon. Not. Roy. Astron. Soc. 363, 



19 



529 (2005) [arXiv:astro-ph/0411170l ; 

K. Takahashi, K. Ichiki, H. Ohno and H. Hanayama, "Magnetic field generation from cosmological perturbations," Phys. 
Rev. Lett. 95, 121301 (2005) [arXiv:astro-ph/0502283] ; 

K. Ichiki, K. Takahashi, H. Ohno, H. Hanayama and N. Sugiyama, "Cosmological Magnetic Field: a fossil of density 



perturbations in the early universe," arXiv:astro-ph/0603631 



K. Ichiki, K. Takahashi, N. Sugiyama, H. Hanayama and H. Ohno, "Magnetic Field Spectrum at Cosmological 
Recombination," |arXiv:astro-ph/0701329 



T. Kobayashi, R. Maartens, T. Shiromizu and K. Takahashi, "Cosmological magnetic fields from nonlinear effects," 
|arXiv:astro-ph/0701596| 

[15] V. Acquaviva, N. Bartolo, S. Matarrese and A. Riotto, "Second-order cosmological perturbations from inflation," Nucl. 



Phys. B 667, 119 (2003) arXiv:astro-ph/0209156 . 
[16] U. Seljak and M. Zaldarriaga, "A Line of Sight Approach to Cosmic Microwave Background Anisotropics, " Astrophys. J. 

469, 437 (1996) [arXiv:astro-ph/9603033] ; |http : //www. cmbf ast . org| 
[17] A. Lewis and A. Challinor; http://camb.info 

[18] W. Hu, D. Scott, N. Sugiyama and M. J. White, "The Effect Of Physical Assumptions On The Calculation Of Microwave 

Background Anisotropics," Phys. Rev. D 52, 5498 (1995) arXiv:astro-ph/9505043] . 
[19] S. Weinberg, "Damping of tensor modes in cosmology," Phys. Rev. D 69, 023503 (2004) [arXiv:astro-ph/0306304] ; 

Y. Watanabe and E. Komatsu, "Improved calculation of the primordial gravitational wave spectrum in the standard 

model," Phys. Rev. D 73, 123515 (2006) [arXiv:astro-ph/0604176l . 
[20] David Spergel and Lyman Page, private communication. 
[21] Latham Boyle, PhD Thesis, Princeton University (2006). 
[22] D. Baumann and K. Takahashi, unpublished. 

[23] T. L. Smith, M. Kamionkowski and A. Cooray, "Direct detection of the inflationary gravitational wave background," 

Phys. Rev. D 73, 023504 (2006) [arXiv:astro-ph /0506422] . 
[24] D. Baumann, K. Ichiki, P. Steinhardt, K, Takahashi, in preparation. 



